Quantum light by atomic arrays in optical resonators 



O 

(N 

< 
(N 



Oh- 



> 
00 

o 

O 



X 



Hessam Habibian, 1 ' 2 Stcfano Zippilli, 1, 2 ' 3 and Giovanna Morigi 1,2 

1 Grup d 'Optica, Departament de Fisica, Universitat Autdnoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain 
2 Theoretische Physik, Universitat des Saarlandes, D-66041 Saarbriicken, Germany 
1 Department of Physics, Technische Universitat Kaiserslautern, D-67663 Kaiserslautern, Germany 

(Dated: January 18, 2013) 

Light scattering by a periodic atomic array is studied when the atoms couple with the mode of a 
high-finesse optical resonator and are driven by a laser. When the von-Laue condition is not satified, 
there is no coherent emission into the cavity mode, and the latter is pumped via inelastic scattering 
processes. We consider this situation and identify conditions for which different non-linear optical 
processes can occur. We show that these processes can be controlled by suitably tuning the strength 
of laser and cavity coupling, the angle between laser and cavity axis, and the array periodicity. We 
characterize the coherence properties of the light when the system can either operate as degenerate 
parametric amplifier or as a source of antibunched-light. Our study permits us to identify the 
individual multi-photon components of the nonlinear optical response of the atomic array and the 
corresponding parameter regimes, thereby in principle allowing one for controlling the nonlinear 
optical response of the medium. 

PACS numbers: 42.50.Nn, 42.50.Pq, 42.65. Yj 



I. INTRODUCTION 

Resonance fluorescence from a single atom exhibits 
non-classical features [l|, which become evident in the 
correlation functions of the emitted light @]. Non- 
classical properties emerge from the quantum nature of 
the scatterer, such as the discrete spectrum of the elec- 
tronic bound states of the scattering atom. They can be 
enhanced or suppressed by several scatterers forming a 
regular array [J-Q . In this case, at the solid angles which 
satisfy the von-Laue condition Q , the light in the far field 
is in a squeezed coherent state, while for a large number 
of atoms it can exhibit vacuum squeezing at scattering 
angles, for which the elastic component of the scattered 
light is suppressed [j§. 

When the atoms of the array are strongly coupled with 
the mode of a high-finesse resonator, emission into the 
cavity mode is in general expected to be enhanced. The 
properties of the light at the cavity output will depend on 
the phase-matching conditions, determined by the angle 
between laser and cavity wave vector and by the peri- 
odicity of the atomic array. The coherence properties 
of the light at the cavity output may however be signif- 
icantly different from the ones predicted in free space. 
An interesting example is found when the geometry of 
the setup is such that the atoms coherently scatter light 
into the cavity mode. In this case the intracavity-field 
intensity becomes independent of the number of atoms 
N as N increases, while inelastic scattering is suppressed 
over the whole solid angle in leading order in \/N Q. 
These dynamics have been confirmed by experimental 
observations (M [t| , and clearly differ from the behaviour 
in free space Q. 

When the geometry of the setup is such that the von- 
Laue condition is not satisfied, photons can only be in- 
clastically scattered into the cavity mode. The smaller 
system size for which coherent scattering is suppressed 



is found for two atoms inside the resonator. The prop- 
erties of the light at the cavity output for this specific 
case have been studied in Refs. [13: El- To the best of 
our knowledge, however, the scaling of the dynamics with 
the number of atoms N is still largely unexplored in this 
regime. 

In this article we characterize the coherence properties 
of the light at the cavity output when the light is scat- 
tered from a laser into the resonator by an array of atoms 
and the geometry of the system is such that coherent 
scattering is suppressed. For the phase-matching condi- 
tions, at which in free space the light is in a squeezed- 
vacuum state 0, we find that inside a resonator and 
at large N the system behaves as an optical paramet- 
ric oscillator, which in certain regimes can operate above 
threshold [l2||- For a small number of atoms N, on the 
contrary, the medium can act as a source of antibunched 
light. In this case it can either behave as single-photon 
or, for the saturation parameters here considered, two- 
photon "gateway" [l3| . The latter behaviour is found for 
a specific phase-matching condition. We identify the pa- 
rameter regimes which allow one to control the specific 
nonlinear optical response of the medium. 

This article is organized as follows. In Scc.|TI]thc theo- 
retical model is introduced and the basic approximations 
are discussed, which lead to the derivation of an effective 
Hamiltonian for the atomic and cavity excitations. In 
Sec. IIIII we analyze the light at the cavity output under 
the condition that there is no coherent scattering into 
the cavity. The conclusions arc drawn and outlooks are 
discussed in Sec. [TV] 



II. THEORETICAL MODEL 

The physical system is composed by N identical atoms 
which are regularly distributed along the z axis. They 
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FIG. 1: (Color online) (a) An array of atoms, with inter- 
particle distance d, is confined along the axis of a standing- 
wave optical cavity at frequency u c and is transversally driven 
by a laser, whose wave vector forms the angle O with the 
cavity axis. The atomic internal transition and the relevant 
frequency scales are given in (b), with |1) and |2) ground 
and excited state of an optical transition with frequency coo 
and natural linewidth 7. The frequencies uj z = u)a — 0J P and 
S c — ui c — ui p denote the detunings between the laser fre- 
quency uip and the atomic and cavity frequency, respectively. 
The other parameters are the laser Rabi frequency f2, the 
atom-cavity coupling strength g, and the decay rate k of the 
optical cavity. 



are located at the positions Zj = jd w here j = 1, . . . , N 
and d is the interparticle distance [14| . An optical dipole 
transition of the confined atoms interacts with the mode 
of a standing wave cavity, whose wave vector k is parallel 
to the atomic array, as illustrated in Fig. Q] Moreover, 
the atoms are transversally driven by a laser and scatter 
photons into the cavity mode. We denote by lu p the fre- 
quency of the laser mode, whose polarization is assumed 
to be linear and parallel to the polarization of the cavity 
mode, and a and a* the annihilation and creation op- 
erators, respectively, of a cavity photon at frequency uj c 
(with [a,al] = 1). Cavity and laser mode couple to the 
atomic dipolar transition at frequency loq with ground 



and excited states |1) and \2). The Hamiltonian govern- 
ing the coherent dynamics of cavity mode and the atoms 
is given by 



N 



TL =hw c a) a + Hluq Sj 



3=1 



JV 



fig cos (kzj + ip)(Sja + a* 1 Sj 



0=1 

N 

+ iM^(Sje- i ^*e i(fc ^ cose - < ^ ) -H.c), (1) 
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where f2 is the strength of the coupling between laser 
and atomic transition, and g the cavity vacuum Rabi 
frequency. The operators Sj = |l)j(2| and Sj indicate 
the lowering and raising operators for the atom at the 



ii)jai 



is the z com- 



position Zj, and S? = ^(\2}j(2\ 

ponent of the pseudo-spin operator. In Eq. ([T]) we have 
introduced the angle ip, which is the phase offset of the 
standing wave at the atomic positions, the phase of the 
laser (pL , and the angle between the laser and the cav- 
ity wave vector. 

We remark that in the present study we neglect the 
atomic motion, and consider that the size of the atomic 
wavepacket is much smaller than the laser wavelength 
and interparticle distance. We refer the reader to [Hj 
for a quantitative study on the effects of the mechanical 
motion on the nonlinear optical processes in this kind of 
system. 

In the rest of this section we will introduce and discuss 
the approximations, which allow us to solve the dynam- 
ics and determine the properties of the cavity field. For 
simplicity we also set k = k p : The difference between the 
laser and cavity wave numbers can in fact be neglected 
for the purpose of this study. 



A. Weak excitation limit 

We resort to the Holstein-Primakoff representation for 
the spin operator [15 1 



Sj = (1 - bjbj^bj, 



(2) 



where bj (bj) is the bosonic operator annihilating (creat- 
ing) an excitation of the atom at Zj, such that [bj, bj,] = 
Sjji . In the limit in which the atomic dipoles are driven 
below saturation, we treat saturation effects in the lowest 
non-vanishing order of a perturbative expansion, whose 
small parameter is the total excited-state population of 
the atoms, denoted by A/tot- We denote the detuning of 
the laser from the atomic transition by 



^0 ~ w p 



(3) 
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and by 7 the natural lincwidth. In the low saturation 
limit, \u) z + vy/2\ > y/NCl, then Af tot < N and we 
can expand the operators on the right-hand side of the 
equations @ in second order in the small parameter 
(bjbj) 1, obtaining 



(4) 
(5) 



For N 3> 1 the dynamics is expected to be irrelevantly 
affected by the assumptions on the boundaries. There- 
fore, we take periodic boundary conditions on the lattice, 
such that zn+i = z\. The atomic excitations are studied 
in the Fourier transformed variable q, quasimomentum 
of the lattice, which is defined in the Brillouin zone (BZ) 
q € (—Go/2, Go/2] with Go = 2n/d the primitive recip- 
rocal lattice vector. Correspondingly, we introduce the 
operators b q and b q , defined as 
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(6) 



(7) 



which annihilate and create, respectively, an excitation 
of the spin wave at quasimomentum q and fulfilling the 
commutation relation [b q , b ql ] = 5 q . q > . After rewriting the 
Hamiltonian in Eq. (fTJ) in terms of spin- wave operators, 
we find 



H 



Nhjj r 



(8) 



where the first term on the Right-Hand Side (RHS) is a 
constant and will be discarded from now on, while 



i(bj p t+<p L ) 



(9) 



is the linear term describing the coupling with the laser. 
Term 



Tl^ =fko c a^a + Huq 
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H.c. 



(10) 



determines the system dynamics when the linear pump is 
set to zero and the dipoles are approximated by harmonic 
oscillators (analog of the classical model of the elastically 
bound electron), while 
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(ii) 



accounts for the lowest-order corrections due to satura- 
tion. In Eqs. (fTU|) and (TTT|) we have denoted by ±Q and 
Q' the quasimomenta of the spin waves which couple to 
the cavity and laser mode, respectively, and which fulfill 
the phase matching conditions 



Q = k + G, 

Q' = k cos 9 + G' , 



(12) 
(13) 

with reciprocal vectors G,G' such that Q, Q' £ BZ. The 
atoms scatter coherently into the cavity mode when the 
von-Laue condition is satisfied, namely one of the two 
relations is fulfilled: 



2fcsin 2 (9/2) = nG 
2/ccos 2 (6/2) = n'G 



(14) 
(15) 

with n, n' integer numbers. In free space, the von-Laue 
condition corresponds to Eq. (|14[) : for these angles one 
finds squeezed-coherent states in the far field When 
the scattered mode for which the von-Laue condition 
is fulfilled corresponds to a cavity mode, superradiant 
scattering enhances this behaviour, until the number of 
atoms N is sufficiently large such that the cooperativity 
exceeds unity. In this limit one observes saturation of 
the intracavity-ficld intensity, which reaches an asymp- 
totic value whose amplitude is independent of N as N 
is further increased. In the limit N ^> 1 the light at 
the cavity output is in a coherent state, while inelastic 
scattering is suppressed at leading order in 1/N @. 

When the von-Laue condition is not satisfied, classi- 
cal mechanics predicts that there is no scattering into 
the cavity mode. These modes of the electromagnetic 
fields are solely populated by inelastic scattering pro- 
cesses. Moreover, in free space, when 



2fcsin 2 (e/2) = (2n + l)G /2 : 



(16) 



then the inelastically scattered light is in a vacuum- 
squeezed state [||. Inside a standing- wave resonator, on 
the other hand, the mode is in a vacuum-squeezed state 
provided that cither Eq. (I16|) or an additional relation, 

2fccos 2 (6/2) = (2n+ l)G /2, (17) 

is satisfied. 

In the following we will study the field at the cavity 
output as determined by the dynamics of Hamiltonian ([8]) 
when Q' ^ ±Q, namely, when the scattering processes 
which pump the cavity are solely inelastic. We remark 
that throughout this treatment we do not make specific 
assumptions about the ratio between the array periodic- 
ity d and the light wavelength A (and therefore also con- 
sider the situation in which A 7^ 2d. This situation has 
been experimentally realized for instance in Refs. [H, liH - 

El). 



B. Linear response: Polaritonic modes 

We first solve the dynamics governed by Hamiltonian 
'HP^ in Eq. (JTUJ). In the diagonal form the quadratic part 
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can be rewritten as 

2 



n^=j2^nh+ E ^ob\b q , (is) 



q^Q s ,qeBZ 



where Q s labels the spin wave which couples with the 
cavity mode, such that 



if Q = 0,G /2, 



(19) 



bQs = S9l - otherwise . (20) 



The resulting polaritonic eigenmodes are 

71 = —a cos X + bg^ sin X , 

72 = a sin X + 6q s cos X , 

with respective eigenfrequencies 

oj 1:2 = - (oj c + cj T <M , 



and 



tan X = gv N / (ujq — lj\ ) , 



(21) 
(22) 



(23) 
(24) 



(25) 



which defines the mixing angle X. The parameter g is 
proportional to the coupling strength. In particular, g — 
gcosip when Q = 0,Go/2 and the cavity mode couples 
with the spin wave &q s = 6g, otherwise g = g/V% and 
the spin wave is given in Eq. ([2"U]) . 

Hamiltonian ^ describes the coupling of the pump 
with the spin wave Q' . When Q' ^ ±Q, photons are 
pumped into the cavity via inelastic processes, which in 
our model are accounted for by the Hamiltonian term in 
Eq. (fTTj) . On the other had, when the dynamics is con- 
sidered up to the quadratic term (hence, inelastic pro- 
cesses are neglected), only the mode Q' is pumped and 
the Heisenberg equation of motion for &q< reads 



-be 



+ nVxc-^ + ViKUt) , (26) 



that has been written in the reference frame rotating at 
the laser frequency uj p . Here, 7 is the spontaneous de- 
cay rate and b qj i n (t) is the corresponding Langevin force 
operator, such that (b q s n (t)) = and (b q ,i n (t)b q in (t')) = 
5(t — t') [l2j ■ The general solution reduces, in the limit 
in which \u> z \ ^> 7, to the form 



'Q' 



(27) 



which is consistent with the expansion to lowest order 
in Eq. provided that fl 2 N <C u 2 - In the reference 
frame rotating at the laser frequency the explicit fre- 
quency dependence of the Hamiltonian terms is dropped, 
and ui\ — > lui — uj p . uj 2 ^2 — ^p, ^0 ~^ w z , and 

UJ r — > 0J r — LO n = S r - 



C. Effective Hamiltonian 

Under the assumptions discussed so far, we derive from 
Hamiltonian ^ an effective Hamiltonian for the polari- 
ton 71 . The effective Hamiltonian is obtained by adiabat- 
ically eliminating the coupling with the other polaritons, 
according to the procedure sketched in Appendix A, and 
reads 



H 



eff 



=Mwi7|7i 

1 ft I \2 2id> L 1 * 2 -2icb L \ , fc ft 

+ 2 \ a n e + a 7i e ) + frxjiihui 



where ll£ 



2 ^ 2 , 

ou)\ — u)\ — uj v H I b 



2 , m S e 



(28) 



(29) 



--(S 2 + ^SC) (<W + ^ G/2 ) (30) 



1+C k^G/2 



9 S 3 C 



(31) 



with S = sin A and C 



S 2 C C^ G/2 , (32) 



cosX. The terms £ fe ^ G / 2 do 



not vanish when k ^ G/2, and their explicit form is 



,^G/2 _ 



r k^G/2 



c 



k^G/2 



\ + ^Q,±Go/4 COS (4(f)j (1 - S KG/2 ), 

\ {h',Q+G/2e- 2ltp + S Q ,^ Q+G/2 e 2 ^) 

X(l - 4,G/2), 

^= (s Q ',3 Q e- 3 ^ + 8 Q ',-iQ^) (1 ~ Kg/2 



The coefficients have been evaluated under the require- 
ment Q' ^ ±Q. 

We now comment on the condition £l 2 N <^ uj 2 , on 
which the validity of Eq. (|27|) is based. When this is 
not fulfilled, such that |(&q')I ~ 1) Eq. (|2^|) must con- 
tain further non-anharmonic terms from the expansion 
of Eq. ([3]) and which account for the saturation effects in 
bqi . Since this spin mode is weakly coupled to the other 
modes, which are initially empty, we expect that the po- 
laritons 71 and 72 will remain weakly populated and the 
structure of their effective Hamiltonian will qualitatively 
not change. 

We also note that in the resonant case, when ui z = 0, 
the form of Hamiltonian in Eq. (|28p remains unchanged, 
while in the coefficients <5wi, a, v, x the- following sub- 
stitution ui z — > (lo 2 + 7 2 /4)/w z must be performed . A 
first consequence is that a = 0, which implies that there 
are no processes in this order for which polaritons are 
created (annihilated) in pairs. A further consequence is 
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that spontaneous decay plays a prominent role in the dy- 
namics. We refer the reader to Sec. Ill El for a discussion 
of the related dissipative effects. 



D. Discussion 

Let us now discuss the individual terms on the RHS of 
Eq. (f28|). For this purpose it is useful to consider multi- 
level schemes, which allow one to illustrate the relevant 
nonlinear processes. The multilevel schemes are depicted 
in Fig. [3] state \n) is the polariton number state with n 
excitations. The blue arrows indicate transitions which 
are coupled by the laser, for which the polariton state 
is not changed; The red arrows denote transitions which 
are coupled by the cavity field, for which the polariton 
state is modified by one excitation. 

Using this level scheme, one can explain the dynamical 
Stark shift 5u>i of the polariton frequency in Eq. (|2U)) as 
due to higher-order scattering processes, in which laser- 
and cavity- induced transition creates and then annihi- 
lates, in inverse sequential order, a polariton. 

The second term on the RHS has coupling strength 
given in Eq. p0[) . it generates squeezing of the polari- 
ton and does not vanish provided that Q' = G/2 or 
Q' = ±Q + G/2. The latter condition is equivalent to 
the free-space condition (|f 6[) . while the first arises from 
the fact that the cavity mode couples with the symmet- 
ric superposition 6q s in Eq. (TT9)) . The corresponding 
phase-matched scattering event is a four-photon process, 
in which two laser photons are absorbed (emitted) and 
two polaritonic quanta are created (annihilated). For 
Q' = G/2 and Q' ^ ±Q the polaritons are created 
in pairs with quasimomcntum Q and — Q (the relation 
Q = G/2 corresponds to 6_q = bq). This specific term 
is also present when the geometry of the setup is such 
that von-Laue condition is fulfilled, and at this order is 
responsible for the squeezing present in the light at the 
cavity output. 

The Kerr-nonlinearity (third term on the RHS) gives 
rise to an effective interaction between the polaritons and 
emerges from processes in which polaritons are absorbed 
and emitted in pairs. It is depicted in Fig. 2(b) for a 
generic case. This term is directly proportional to the 



cavity coupling strength and inversely proportional to 
y/W. In this order, it is the term that gives rise to anti- 
bunching. 

The last term on the RHS, finally, is a nonlinear pump 
of the polariton mode, whose strength depends on the 
number of polaritonic excitations. It is found when the 
phase-matching condition Q' = ±3Q + G is satisfied, 
which is equivalent to the relation cos (9 = (3 + nX/d) 
when laser and optical resonator have the same wave- 
length, as in the case here considered. This relation can 
be fulfilled for n ^ and specific ratios X/d. This term 
vanishes over the vacuum state, and it pumps a polari- 
ton at a time with strength proportional to the number 
of polariton excitations. 



(a) For Q' = G/2, Q' = 


±Q + G/2 




\h - 2} 






\h-l) 


/ |r* — 2> 









|n) 


l«- 2) 



(a) 




FIG. 2: (Color online) Schematic diagram of transitions which 
fulfill the phase-matching condition till fourth order. State |n) 
denote the number state of the polariton mode 71 . The blue 
(light gray) arrows denote the laser-induced couplings and the 
red (dark gray) arrows denote the creation (annihilation) of 
polaritons due to the coupling with the cavity field. See text 
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In general, photons into the cavity mode are pumped 
provided that either (i) Q' = ±Q or (ii) (for Q' ^ ±Q) 
one of the two conditions are satisfied: Q' = G/2 or 
Q' = ±Q + G/2. We note that the strength of the Rabi 
frequency and of the cavity Rabi coupling may allow 
one to tune the relative weight of the various terms in 
Hamiltonian (|2"5)) . Their ratio scales differently with the 
number of atoms in different regimes, which we will dis- 
cuss below. Moreover, the interparticlc distance of the 
atomic array constitutes an additional control parame- 
ter over the nonlinear optical response of the medium. 
Further phase-matching conditions are found when con- 
sidering higher-order terms in the expansion of the spin 
operators in harmonic-oscillator operators from Eq. ([2]). 
Their role in the dynamics will be relevant, as long as 
they compete with the dissipative rates, here constituted 
by the cavity loss rate and spontaneous emission. 

E. Cavity input-output formalism 

We consider the full system dynamics, including the 
atomic spontaneous emission and the cavity quantum 
noise due to the coupling to the external modes of the 
electromagnetic field via the finite transmittivity of the 
cavity mirrors. Denoting by k and 7 the linewidth of the 
cavity mode and of the atomic transition, respectively, 
the Heisenberg-Langevin equations for the operator a and 
b q read 

a = — [a, H] — na + V2~Ka in (t) , (33) 

b\ = i [b q , H] - \b q + V7VinW , (34) 

where a; n , b q ^ n are the Langevin operators, which are 
decorrelated one from the other and fulfill the relations 
(a in (t)) = and (a in (t)al n (t')) = 6{t - f). Here, the 
average (•) is taken over the density matrix at time t = 
of the system composed by the atomic spins and by the 
electromagnetic field. The output field a ou t at the cavity 
mirror is given by the relation 

a ou t + a in = y/2m(t) . (35) 

Let us now consider the scattering processes occurring 
in the system. They can be classified into three types: 

(i) a laser photon can be scattered into the modes of the 
external electromagnetic field (emf) by the atoms, with- 
out the resonator being pumped in an intermediate time; 

(ii) a laser photon can be scattered into the cavity mode 
by the atom and then dissipated by cavity decay; (iii) a 
laser photon can be scattered into the cavity mode by the 
atom, then been reabsorbed and emitted into the modes 
of the external emf. Processes of kind (i) include elas- 
tic scattering. They can be the fastest processes, but do 
not affect the properties of the light at the cavity out- 
put. Processes of kind (ii) are the ones which outcouple 
the intracavity field, but need to be sufficiently slow in 



order to allow for the build-up of the intracavity field. 
Processes of kind (iii) are detrimental for the nonlinear 
optical dynamics we intend to observe, as they introduce 
additional dissipation (see for instance j2(J [21| for an ex- 
tensive discussion and [l(| for a system like the one here 
considered but composed by two atoms). 

Processes (iii), i.e., reabsorption of cavity photons fol- 
lowed by spontaneous emission, can be neglected assum- 
ing that the laser and cavity mode are far-off resonance 
from the atomic transition. In this limit, the cavity is 
pumped by coherent Raman scattering processes and an 
effective Heisenberg-Langevin equation for the polariton 
71 can be derived assuming that its effective linewidth 
ki = kcos 2 X + (7/2) sin 2 X fulfilling the inequality 
k\ <C Suj (that corresponds to the condition for which 
the vacuum Rabi splitting is visible in the spectrum of 
transmission (l7l. Ha 122^24) ). We find 

7i = 7fc[7ii%eff] - «i7i + V2~KCa in (t) + v /756 9i in(t)C36) 

which determines the dissipative dynamics of the polari- 
ton. The field at the cavity output is determined using 
the solution of the Heisenberg Langevin equation (f3"6]l 
with Eqs. ([2"I ]) -([2"2 ]) in Eq. ([55]). In some calculations, 
when appropriate we solved the corresponding master 
equation for the density matrix of the polaritonic modes 
71 and 72. 

Some remarks arc in order at this point. Nonlinear- 
optical effects in an atomic ensemble, which is reso- 
nantly pumped by laser fields, have been studied for in- 
stance [25j, where the nonlincarity is at the single atom 
level and is generated by appropriately driving a four- 
level atomic transitions |26[. 

It is important to note, moreover, that Eq. (|36|) is valid 
as long as the loss mechanisms occur on a rate which 
is of the same order, if not smaller, than the inelastic 
processes. This leads to the requirement that the atom- 
cavity system be in the strong-coupling regime. 

III. RESULTS 

We now study the properties of the light at the cav- 
ity output as a function of various parameters, assuming 
that Q 1 = G/2 and that the relations Q 1 ^ ±Q and 
Q' 7^ ±3Q + G hold. Under these conditions the effective 
Hamiltonian in Eq. (|28|) contains solely the squeezing and 
the Kerr-nonlinearity terms, while v = 0. Moreover, we 
assume the condition k± ~ k > 7. 

The possible regimes which may be encountered can 
be classified according to whether the ratio 

e = \a/x\ 

is larger or smaller than unity. In the first case the 
medium response is essentially the one of a parametric 
amplifier. In the second case the Kerr non-linearity dom- 
inates, and polaritons can only be pumped in pairs pro- 
vided that the emission of two polaritons is a resonant 
process. 
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Let us now focus on the regime in which the system 
acts as a parametric amplifier, namely, e » 1. In this 
case one finds that the number of photons at the cavity 
output at time t is 



with 



1 



(7i f 7i)t = 7, 



(a>tut a out)t - 2KC 2 (-f\j 1 ) t 



Il4 sinh 2 (crf) 



2k{ -a 2 

~ 2Klt * ki cosh(2crf) + a sinh(2crf) 

1 - Ki 5 5 



(37) 



Depending on whether a > K\ or a. < Ki, one finds that 
the dynamics of the intracavity polariton corresponds to 
a parametric oscillator above or below threshold, respec- 
tively. In the following we focus on the case below thresh- 
old and evaluate the spectrum of squeezing. We first 



observe that the quadrature x 



(0) 



minimum variance for 9 — it /A and reads 



7ie" ie + jle' e has 



(A. 



«1 



(38) 



where the subscript st refers to the expectation value 
taken over the steady-state density matrix. The squeez- 
ing spectrum of the maximally squeezed quadrature is 



Sov.t(w) =1 + 



+ 00 



(:xlil(t + T),x£>(t) :) st e 



AnC 2 \a\ 



{ki + \a\) 2 +oj 2 



(39) 
(40) 



where (: :) st indicates the expectation value for the 
normally-ordered operators over the steady state, with 



,t 



and (A,B) st = (AB) st - (A) st (B) st . 

We now discuss the parameter regime in which these 
dynamics can be encountered. The relation s 3> 1 is 
found provided that SI ^> g. When ^> Qy/N, in this 
limit |a| ~ Q?g 2 N/u%, and squeezing can be observed 
only for very small values of k. Far less demanding pa- 
rameter regimes can be accessed when relaxing the con- 
dition on the laser Rabi frequency, and assuming that 
fly/N ~ |cj z |. In this case squeezing in the light at the 
cavity output can be found provided that 11 3> k when 
gVN~\w z \ M- 

Figures |3ja) and (b) display the spectrum of squeez- 
ing when the system operates as a parametric amplifier 
below threshold. Here, one observes that squeezing in- 
creases with N. Comparison between Fig. [3Ja) and [3] 
(b) shows that squeezing increases also as the single- 
atom cooperativity increases (provided the correspond- 
ing phase-matching conditions arc satisfied and the laser 
Rabi frequency 3> <?). These results agree and extend 
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0.8 
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0.4- 



0.2- 
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FIG. 3: (Color online) Squeezing spectrum for the maximum 
squeezed quadrature when k = G/2, Q' — G/2 and Q' 7^ Q 
for tp, 4>l = 0. The parameters are = 200k, lo z — 10 3 k and 
N = 10, 50, 100 atoms (from top to bottom) for (a) g — 4k 
and (b) g = 10k. The detuning S c is chosen such that Suji 
is zero. The curves are evaluated from Eq. (|39p by numeri- 
cally calculating the density matrix of the polariton field for 
a dissipative dynamics, whose coherent term is governed by 
the effective Hamiltonian in Eq. (|28|) . The value g = 4k is 
consistent with the experimental data of Ref. [29I ]. 



the findings in Ref. [Toj , which were obtained for an array 
consisting of 2 atoms. 

Let us now focus on the regime when e <C 1. Here, 
the polaritons may be only emitted in pairs into the res- 
onator. In order to characterize the occurrence of these 
dynamics we evaluate the second-order correlation func- 
tion at zero time delay in the cavity output defined by 

m 



5 (2) (0) 



( a out a out)st 
( a out a out)?t 



(41) 



Function g^(0) quantifies the probability to measure 
two photon at the cavity output at the same time. There- 
fore, subpossonian (supcrpossonian) statistics are here 
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connected to the value of g^ 2 \0) smaller (larger) than 
one, while for a coherent state g^ 2 '(0) = 1. 

Subpossonian photon statistics at the cavity output 
can be found as a result of the dynamics of Eq. (|28p. 
Here, for phase-matching conditions leading to v = and 
a/0, polaritons can only be created in pairs. When the 
Kerr-nonlinearity is sufficiently large, however, the con- 
dition can be reached in which only two polaritons can be 
emitted into the cavity, while emission of a larger number 
is suppressed because of the blockade due to the Kerr- 
term. This is reminiscent of the two-photon gateway re- 
alized in Ref. [l3j . where injection of two photons inside 
a cavity, pumped by a laser, was realized by exploiting 
the anharmonic properties of the spectrum of a cavity 
mode strongly coupled to an atom. In the case analysed 
in this paper, the anharmonicity arises from collective 
scattering by the atomic array, when this is transversally 
driven by a laser. Moreover, we note that the observa- 
tion of these dynamics requires Q^/N <C \u x \,g\N and 
\a\ > k, which reduces to the condition Q, 2 /uj z > k when 

Figure 0Ja) displays <r 2 '(0) as a function of the pump 
frequency ui p for the phase matching conditions giving 
v = and a^O. Function <?' 2 '(0) is evaluated by numer- 
ically integrating the master equation with cavity decay, 
where the coherent dynamics is governed by an effective 
Hamiltonian which accounts for the effect of both polari- 
ton modes and is reported in Eq. (|A1[) in the Appendix. 
Antibunching is here observed over an interval of values 
of u) p , about which the cavity mode occupation has a 
maximum (blue curve in Fig. E]Jb)). The maximum cor- 
responds to the value of uj p for which the emission of two 
polaritons 72 is resonant. Note that the spin-wave exci- 
tation, red curve in Fig. BJb), i s still sufficiently small to 
justify the perturbative expansion at the basis of our the- 
oretical model. Figure 4(c) displays the amplitudes |%|, 
determining the strength of the Kerr-nonlinearity, and 
|a|, scaling the squeezing dynamics, in units of \8uj\ \ and 
as a function of uj p . One observes that for the chosen pa- 
rameters \x\ > \a\. Maximum antibunching is here found 
when the cavity mean photon number is maximum. 

It is important to notice that emission of polaritons in 
pairs is possible when the collective dipolc of the atomic 
array is driven. For fixed values of f2 and g, we expect 
that this effect is washed away as TV is increased: this be- 
haviour is expected from the scaling of the ratio e with TV. 
Taking k = G/2 and mixing angles X -C 1, for instance, 
one finds e ~ a/TV, indicating that the strength of the 
Kerr nonlincarity decreases relative to the coupling a as 
TV grows. This is also consistent with the results reported 
in Fig. [3] In this context, the expected dynamics is rem- 
iniscent of the transition from antibunching to bunching 
observed as a function of the number of atoms in atomic 
ensembles coupled with cavity QED setups (28j |. 

For the results here presented we have assumed the 
spontaneous emission rate to be smaller than k. In gen- 
eral, the predicted nonlinear effects can be observed in 
cavities with a large single-atom cooperativity and in the 




134 136 

(UJQ - UJ p )/K 

(a) 




(cj - u p )Jk 
(b) 




134 136 

[UO —Up)/K 

(C) 



FIG. 4: (Color online) (a) Second order correlation function 
at zero time delay <7' 2 '(0) versus oj p (in units of k) when the 
cavity is solely pumped by inelastic processes (here, k — G/2, 
Q' Q,3Q and Q' = G'/2 for <p,<f> L = 0). The correlation 
function is evaluated numerically solving the master equation 
for the polaritons in presence of cavity decay, with the co- 
herent dynamics given by Hamiltonian (|A1|) (solid line) and 
by Hamiltonian (|28[) (dashed line), (b) Corresponding av- 
erage number of intracavity photons (a^a) (blue upper line) 
and spin wave occupation (6q&q) (red lower line), (c) Ratios 
(blue upper line) and |a/<Jwi| (red lower line) versus 



The parameters are g = 



Q = 30k, 



5 C = 70k, 



and TV = 2 atoms. At the minimum of g^ 2 \0), u z ~ 137k and 
5 C ~ 67k. 
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so-called good cavity regime [22[. The required param- 
eter regimes for obser ving squeezing have been realized 
in recent experiments [20|. The parameters required in 
order to observe a two-photon gateway are rather de- 
manding for the regime in which the atoms are driven 
well below saturation. Nevertheless, a reliable quantita- 
tive prediction with an arbitrary number of atoms would 
require a numerical treatment going beyond the Holstcin- 
PrimakofF expansion here employed. 



with radiation [lflj . It is important to consider, that when 
the mechanical effects of the scattered light on the atoms 
is taken into account, conditions could be found where 
selforganized atomic patterns are observed [HI, |H, HH- 
[39j . which are sustained by nonclassical light. The anal- 
ysis here presented sets the basis for studies towards this 
direction. 
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IV. CONCLUSION 

An array of two-level atoms coupling with the mode 
of a high-finesse resonator and driven transversally by a 
laser can operate as controllable nonlinear medium. The 
different orders of the nonlinear responses correspond to 
different nonlinear processes exciting collective modes of 
the array. Depending on the phase-matching condition 
and on the strength of the driving laser field a nonlinear 
process can prevail over others, determining the dom- 
inant nonlinear response. These dynamics emerge from 
the backaction of the cavity mode on the scatterers prop- 
erties, and are enhanced for large single-atom coopcra- 
tivities. We have focussed on the situation in which the 
scattering into the resonator is inelastic, and found that 
at lowest order in the saturation parameter the light at 
the cavity output can be either squeezed or antibunched. 
In the latter case, it can either operate as single-photon or 
two-photon gateway, depending on the phase-matching 
conditions. Our analysis permits one to identify the pa- 
rameter regimes, in which a nonlinear-optical behaviour 
can prevail over others, thereby controlling the medium 
response. An interesting outlook is whether the consid- 
ered effects can be used in order to develop turnstile de- 
vices, like the one realised in [30j . for quantum networks. 

In view of recent experiments couplin g ul tracold atoms 
with optical resonators @, i, G3, H HIM EMI , these 
findings show that the coherence properties at the cav- 
ity output can be used for monitoring the spatial atomic 
distribution inside the resonator. Another related ques- 
tion is how the properties of the emitted light depend 
on whether the atomic distribution is bi- or multiperi- 
odic (34[. In this case, depending on the characteristic 
reciprocal wave vectors one expects a different nonlin- 
ear response at different pump frequency and possibly 
also wave mixing. When the interparticle distance is 
uniformly distributed, then coherent scattering will be 
suppressed. Nevertheless, the atoms will pump inclasti- 
cally photons into the cavity mode. While in free space 
the resonance fluorescence is expected to be the incoher- 
ent sum of the resonance fluorescence from each atom, 
inside a resonator one must consider the backaction due 
to the strong coupling with the common cavity mode. 

Finally, in this article we neglect the atomic kinetic en- 
ergy, assuming that the spatial fluctuations of the atomic 
center of mass at the potential minima are much smaller 
than the typical length scales determining the coupling 
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Appendix A: Derivation of the effective Hamiltonian 

^From the general form of the Hamiltonian in Eq. (|SJ) 
for the case in which there is no coherent scattering 
(Q' 7^ ±Q), in the weak excitation limit one can obtain 
the effective dynamics for the polariton described by 71. 
We focus on the regime in which ily/N <C |w z |. As we 
are interested in the dynamics of the mode 6q s and of 
the cavity mode a, the relevant terms determining their 
dynamics are given in lowest order by 

(7 = 1,2 



(Al) 



with 

w 



ft (J 



4v/A 



{btybty (b. Q e^ + b Q e-^) 8 Q , >G/ 2 



+2bl (b^ + blqe-*) b Q 



blblbne^ 



b-nC~ lv 



+(1 - <W) [26^6Lq6-q + *q,±g /4& W« 

+<*Q',Q+G/2&Q'&Q'&Q e' V 



-(1 



■>k,G/2 



2tf Q bl Q b Q + S Q>±Go/4 b J Q b T Q b 



t it 



+2^- 3 g,g'&q 6 g' 6 -Q + s 3Q,Q' b Q b Q b Q 



t At 



~ [ i7N e ^ L { 2b Q' b Q bQ + 2(1 ~ s ^ b -Q b h b -> 

+SQ>,G/2b* Q bl Q b(y [1 + (1 - <5 fc ,G/ 2 )] 
+ (1 - 4,G/2, ^30,0'6q6q6-Q + £-3Q,Q'&1q&1q&<; 
+h',Q+G/2 b Q b Q b Q' +5Q l _Q +G/ 2bl bl b Q i I 

+H.C. 



By substituting &q< with its mean value in Eq. (|27[) , which 
corresponds to neglect the backaction on the mode Q' due 
to the nonlinear coupling, one obtains closed equations 
of motion for the modes &q s and a (where we thereby 
discard the effect of the nonlinear coupling with the other 
modes, which are initially empty and which gives rise 
to higher order corrections). In this limit the effective 
Hamiltonian (|28|) for the polariton 71 is derived provided 
that the detuning of the laser from the polariton 72 is 
much larger than the strength of the nonlinear coupling 
with polariton 71. 
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